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1. Large QCD in the 't Hooft limit 

Large N gauge theories ^ are qualitatively similar to QCD with three colors and it is long held 
hope to solve it analytically for N = 00. Only planar diagrams contribute in this limit and fermions 
in the fundamental representation are naturally quenched as long as the number of flavors is finite. 
Researchers in string theory and gravity also address the problem of large N gauge theories but 
they are a long way from solving it analytically 

We will discuss various physical properties of large N QCD in the 't Hooft limit. We will start 
with well known results m d = 2 and proceed to a discussion of relatively new results in J = 3 
and d = A. The theory will be regulated using the lattice formalism and the lattice bare coupling, 
b = will be held fixed as g ^ and N —>■ °°. We will assume that we only have a finite 
number of fermion flavors and therefore fermions will be naturally quenched as long as we are in 
the confined phase and there is no chemical potential. The continuum limit corresponds to ^ 00. 
All our discussion will be on a periodic lattice at a finite physical volume. Specifically, 

• d = 2: Ijcy = The physical sizes < ly are kept fixed as L^ y and b are taken to 00. 

• d = 3: lx,y,z — \ 

We will use the tadpole improved coupling |j3|], bj = be{b) where e{b) is the average value 
of the plaquette. The physical sizes 4 < h < h kept fixed as Lvj,^ and b are taken to 00. 



— 4: lx,y,z,t — ^x,y,zJ 

Lx,y,z,t and b are taken to 00. 
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e 11 . The physical sizes < ly <l^<lt are kept fixed as 



2. f/( 1) symmetry and continuum reduction 



The lattice gauge action for SU {N) gauge theory on a Li x • • - L^ periodic lattice is given by 

- £ rr[[/^,v(«) + t/;,v(«)] (2.1) 



2 n,ji^v 



U^,y{n) = Uy,{n)U^{n + pL)Ul{n + v)Ul{n). (2.2) 

In addition to the local gauge symmetry, the above action has a global symmetry under which 
the Polyakov loop in the d directions get rotated by a Za? phase factor: 

U^{n) e''^'"'>'l'^U^{n) for n^=L^ and < ?iv < for ^ / V Q<k^<N (2.3) 

Each Z^r becomes a ?/( 1) in the A/^ — > 00 limit. 

If the U{V) symmetry is not broken in a given direction on a Li x •••Lrf lattice at a fixed 
coupling b, then no physical quantity depends on the size of that direction. The proof of the above 
statement is a simple extension of the original Eguchi and Kawai argument to a Li x •••Lrf 
lattice. The continuum limit of the above statement (namely, 00^ b 00^ such that the physical 
size is kept fixed) is referred to as continuum reduction 

If continuum reduction holds in a certain direction, the parallel transporter in that direction can 
be folded using periodic boundary conditions to construct a transporter of arbitrary length. This 
enables one to consider Wilson loops of arbitrary size on a finite box. 
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Continuum reduction has interesting consequences for fermions. First one notes that the gauge 
transformations can be extended from SU (N) to U(N) and still get the same result for fermionic 



gauge invariant quantities. One can then convert ( |2.3[ ) to 



(n) e lJ^{n)\ 0<k^<N (2.4) 

using a U(N) gauge transformation that obeys periodic boundary conditions. The gauge field ac- 
tion is invariant under the above transformation and therefore an observable made out of a single 
fermion (like the quark condensate) cannot depend upon k. But a fermionic observable that is made 
out of one quark and a different anti-quark (like a 71+ meson) will depend upon {k — q) if one quark 
sees a gauge field with k^ and the other quark sees a gauge field with . Therefore, one can have 
continuous momenta in this direction where the discete momentum interval on the lattice are filled 
by —jji^- This is called the quenched momentum prescription [||] for mesons. 

3. Large QCD in two dimensions 

The two U{\) symmetries remain unbroken for all values of b and L^^y [^. Therefore the 
problem can be reduced to a single site on the lattice with U\ and U2 being the two SU {N) degrees 
of freedom. There is no dependence on the box size or ly for any < Ix < ly < °°- Large N QCD 
in J = 2 is always in the confined phase and there is no dependence on the temperature. 

3.1 Gross- Witten transition 

The plaquette operator is P = U1U2UIU2 and its eigenvalues e'^f", j = l,---,N are gauge 
invariant. Consider the eigenvalue distribution, p{6p;b),for—7l<dp<n obtained upon averaging 
over Ui,2 using the Wilson action. Gross and Witten showed that this observable exhibits non- 
analytic behavior as a function of b. 

p(0,;Z;)=fcos|^^-sin2^,Z,>i; |0,|<2sin-i^ 

p{dp;b)= ^{l+lbcosdp), b<^; \dp\<n (3.1) 

The eigenvalue distribution has a gap for b > ^ and it does not have a gap for b < j. The lattice 
theory has a third order phase transition at b = ^ and this transition is referred to as the Gross- 
Witten transition. This transition is a lattice phenomenon since the location of the transition does 
not scale with the lattice size and the continuum theory is only in the phase where the eigenvalue 
distribution of the plaquette operator has a gap. 

3.2 Wilson loops in large 2d QCD 

The plaquette operator is a Wilson loop whose area goes to zero as one goes to the continuum 
limit. For a physical loop, consider a rectangular Wilson loop of size nxm. The Wilson loop oper- 
ator can be obtained by folding on a single site lattice and is given by W{n,m) 
Let ? = ^ be the parameter that characterizes the dimensional area. Consider a continuum Wil- 
son loop of a fixed area by taking b 00, nm —^ 00 while keeping t fixed. Since, TrW{n,m) = 
[TrW{l, 1)]"'" [§, it is easy to show that TrW(f) = g-i. 
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One can proceed further and get analytical expressions for TrW"(f). In this context, it is useful 
to consider the generating function F{z,t) = 5+^^=1 ^'^''t',, ^^^^ where z is a complex variable. Then, 
F{z,t) satisfies 



_ 2F{z,t) + \ 
2F{z,t)-l 

and the expectation value of powers of Wilson loops are [ |T0| , pl| , 12, 13] 



The expectation value of the distribution of the eigenvalues, e'^, ofW is given by 



p{d,t) 



-ReF(e'^,t) 
n 



(3.2) 



(3.3) 



(3.4) 



3.2.1 Critical behavior of Wilson loops 

The above results imply a critical behavior of Wilson loops as a function of area. The expecta- 



tion value of arbitrary powers of Wilson loops, {TrW"{t)), as given by (3.3) are analytic functions 
of t. Yet, the eigenvalue distribution, p{d,t), exhibits a non-analytic behavior as a function of t 
One way to see this is to ask what p{7i,t) is as a function of t. Setting z 

tF{ 



2F{-l,t) =tanh- 



-l,f 



1 in (3.2) results in 

(3.5) 



F{—l,t) = is always a solution but the non-zero solution for F{—l,t) when ? > 4 is favored. 
Therefore, the critical point is f = 4 and the distribution has a gap for f < 4 and does not have a gap 



for f > 4. It is also clear from the above equation that F{—l,t) oc ^/t — 4 as f ^ 4 [14]. 

A non-trivial critical behavior is observed if one stays at the critical point, t = 4, and asks for 
the behavior of p(0,4) close to 6 = n. If we set z = —e'^' and f = 4 in (3.2), one finds that p{y,4) oc 

0, is proportional to 



y 3 and therefore the number of eigenvalues in an arc of length dy, near y 
4 „ ... . . . , , . 3 



Nyi. Therefore, the level spacing is proportional to N^^ 1 14]. 

This physical transition in Wilson loops as a function of area from weak coupling {t < 4) to 
strong coupling {t > 4) is called the Durhuus-Olesen transition. 

3.2.2 Double scaling limit of the Durhuus-Olesen transition 

The critical behavior of the Wilson loops as a function of area results in a universal function in 
the double scaling limit where one takes t ^4 and 6 ^ n. The double scaling limit can be studied 
by considering 



ON{y,b) 




det +e 



(3.6) 



Yl'j=i Uj is a product of n independently and identically distributed SU (N) matrices, Uj = 
e'^^i . Each Hj is traceless and its entries are independently distributed Gaussian random variables. 



W 

JeH 



We let e ^ and « — > 00 such that the area, | = f = ne^, is kept fixed. 

One can use an integral representation over Grassmann fields and a perturbation expansion in 



£ to show that [15] 



ON{y,b) 



^p^iV[lncoshp-*(2p->f]_ 



(3.7) 



4 



A survey of large N continuum phase transitions 



Rajamani Narayanan 



lattice phase 



No gap in 
plaquette 
distribution 



V 



continuum phase 



Plaquette distribution develops a gap around 9= 



& 

Ml 

s 



Durhuus-Olesen transition in Wilson loop 
as a function of area 



I 

c 
o 
U 



3 
3 
C 

C 

o 
U 



a 

n 

c 
o 
U 



b=l/2 



't Hooft Coupling b=l/g N 



Figure 1: Summary of large QCD ind = 2 



The double scaling limit amounts to taking b ^ I and y ^0. Using the appropriate scaling expo- 
nents obtained in the previous section, we define scaled variables, ^ and a by 

1 

^ 1 

b = \ + ^=a. (3.8) 



3Ny ^' V3]V 

An expansion in results in the following Generalized Airy integral as the universal scaling 
function: 

lim ON{y,b) = n^,a) = I due-"'-"'''+^" (3.9) 

N^oo J 

The above equation describes the universal behavior in the double scaling limit. The conjecture 
of dimensional reduction is that the universal function i^{^,a) defined in the double scaling limit 
for large N QCD in J = 2 is also obeyed by large N QCD ind = 3 and d = 4. 

We end this section with a pictorial summary of large A'^ 2d QCD in Fig. [I] 

4. Large QCD in three dimensions 

The U^{1) symmetries are spontaneously broken on finite lattices and Eguchi-Kawai reduction 
does not hold. The continuum theory can exist in several phases labeled as Oc, Ic, 2c and 3c 
corresponding to the number ofU{l) symmetries that are broken [||, [16, 17 1. The theory is confined 
in the Oc phase and deconfined in the Ic phase. 
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4.1 Transition in the plaquette operator 

Like in 2d, the eigenvalue distribution of the plaquette operator undergoes a transition from 
having no gap for small b to having a gap for large b JT^]. Numerical work indicates that this 
transition is either second or third order and occurs at Z> « 0.43 for the Wilson gauge action. Despite 
some similarities with the Gross-Witten transition, there is no evidence that the transition in J = 3 
is in the same universality class as Gross-Witten. Like in 2d, the location of the transition does not 
scale with the lattice size and therefore it is a lattice transition. Like in 2d, the continuum theory is 
always in the phase where the eigenvalue distribution of the plaquette operator has a gap. 

4.2 Setting the scale in 2>d large QCD 

Although there is some progress in obtaining analytical results in 3d QCD, much of the results 
are obtained by numerical means. One example where numerical results confront analytical results 
is the case of string tension. Using the Hamiltonian formalism and a parameterization of the spatial 
gauge potential using a complex SL{N,C) matrix, one can analytically obtain a good approximation 



to the vacuum wave-function and thereby obtain a result for the string tension []18|]. The result for 
the string tension is 



1-1^ 



1 



(4.1) 



871^72 

for all A^. 

Lattice computations with N = 2,3,4,5,6,8 have been performed[p^, ^0| ] by computing the 
correlation functions of Polyakov loops at zero spatial momentum. The Polyakov loops themselves 
are constructed using smeared gauge fields. The numerical results have very small statistical errors 
(typically less than 0.4%) and they deviate from the above analytical result. The numerical result 
for the string tension is consistently smaller than the analytical result for all N and the differences 
are large for small N. The extrapolation of the lattice results to — > oo gives 

VaZ7 = 0.1975 ± 0.0002 - 0.0005; (4.2) 

and this has to be compared with the analytical result, = 0.19947114- • • The first error in 



the numerical estimate is statistical and the second error which is always negative comes from 
performing two different fits (either a single cosh or a double cosh, which takes into account the 
presence of an excited state in the correlation function). The difference between analytical and 
lattice results although statistically significant is still quite small. 

4.3 Deconfinement transition in 3d large A'^ QCD 

Consider large A'^ QCD on a torus at a fixed lattice coupling b. An order parameter suitable 
for studying the phase transitions we are interested in is ^ 

Px,y,z — {Px,y,z) 

2 



p 1 

^x,y,z 2L 



El 



^x,yM = fiUi{n + mi). (4.3) 

m=\ 
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The quantity P^ j.^ takes values in the range [0,0.5] on any gauge field background and one chooses 
the X, y and z directions for each configuration such that < Py < P^. 

Fixing L and N, one finds the b\ (L) such that one of the three U{\) symmetries is broken for 
b > bi{L): P^^y^z = 5 for < b],{L) and P, < \, Py^^ = \ for b > bi(L). One finds that bi{L) is 
independent of for large enough N (N = 47 is usually sufficient). Then one finds that the tadpole 
improved critical coupling, b\i{L), scales with L and li = L/b\i{L) = 5.90(47). This shows that 
there are two phases in the continuum limit and l\ defines a physical size such that 

• All three U{\) symmetries are unbroken for I > l\ and there is no dependence on I in this 
phase (Oc). 

• One of the three U{\) symmetries is broken for I <l\ and the theory depends on the size of 
the broken direction in this phase (Ic). 

The critical size l\ does not depend on Zj,^ since the U{\) symmetries are not broken in the ly or 
directions and continuum reduction holds in those directions. Therefore, the system is in 4 x oo^ jf 
Ix < /i . It is natural to identify the finite direction with that of temperature in the deconfined phase. 
Therefore, Oc to Ic is the deconfinement transition and the deconfinement temperature is 

^ = 7-^=0.86(7) (4.4) 

Since continuum reduction holds in the Oc phase in all three directions, there is no temperature 
dependence in physical quantities in the Oc phase. A latent heat measurement is needed to directly 
estabhsh the order of the phase transition in the large N limit. 

Conventional numerical studies of SU {N) ind = 3 on x Lf with Lf = 3,4,5 and L as high as 
48 indicate the following. Both SU(2) and SU(3) gauge theories exhibit a second order deconfine- 



ment transition [|21|]. The case of SU(4) is marginal and the transition is difficult to establish [|22|]. 
The Z4 spin model has continuously varying exponents. Small Lt indicate the transition is first 
order but larger Lf possibly indicate a second order transition. The transition is clearly first order 



from SU(5) onwards [23|. The large A'^ limit obtained from extrapolating the A/^ = 4,5,6 results [24] 



for the critical temperature is consistent with the critical size for the Oc to Ic transition. 

4.4 Transition to large A'^ QCD in a small box 

Now consider large A'^ QCD in the Ic phase by picking a box of size Lx x L,, x with Lx < 
Ly < and b > b\{Lx). The box size has been chosen such that the U{\) symmetry in the x- 
direction is broken. As b is increased, the U (1) in the Ly direction will break at some b2{Ly,Lx). 
For the special case of Ly = L^, one of the two U{1) will break and the broken direction will be 
called Ly. 

The theory is in the 2c phase for b > b2{Ly,Lx) and the 2c phase exists in the continuum theory 
since b2i{Ly,Lx) scales with Ly. There is a characteristic size associated with the Ic to 2c transition, 
namely, hih), obtained by taking the limit of Ly/b2i{Ly,Lx) as Ly goes to infinity while keeping 
Lx/Ly fixed. This critical size does not depend on since > hilx) and the U{1) symmetry in that 
direction is not broken. Therefore, the system is in 4 x lyX 00 while in the 2c phase. It is natural to 
associate the two finite directions with a small periodic box and the infinite direction as time. One 
cannot address confinement in the 2c phase since only one direction has an infinite extent. 
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1 




Figure 2: Phase diagram in the {lx,ly) plane for Ix < ly < h 



The system goes into the 2c phase only if 4 < ly <l2 = 0.65 (9)/i . For < Zji < h, the system 
goes into the 2c phase at some hih) with hih) = h- Numerical extrapolation of hih) indicates 
that /2(0) > 0. The full picture is shown in Fig. ^. Note that the critical line connecting the Ic 
and 2c phase is such that one can start in the 2c phase with = ly < h and go into the Ic phase 
by keeping ly fixed and reducing 4- Finally, rotational symmetry is present in the two broken 
directions if 4 = ly < I2. 

4.5 Large A'^ QCD in a small box at high temperature 

Large N QCD on a x Ly x box with b > b2{Lx,Ly) undergoes a phase transition at 
b3{Lx,Ly,L^) beyond which all three U{1) symmetries are broken. The system is in the 3c phase 
for b > bj{Lx,Ly,L^) and corresponds to large A'^ QCD in a small box at high temperature. 

It should be possible to do perturbation theory deep in the 3c phase but one has to account 
for the zero modes of the gauge fields on the torus. There are no zero modes to deal with if one 
considers the theory on x S^. For a small radius of S^, one can show using perturbation theory 
that the U{1) symmetry associated with is broken when the radius of gets smaller than a 
certain size [p5|]. 

Numerical computations show that bi,i{L,L,L) scales properly with L and therefore the con- 
tinuum theory can also exist in the 3c phase. h{lx,ly) is the characteristic size associated with the 
2c to 3c phase transition and h{h,h) = h= 0.36(5)/i. 

4.6 Wilson loop operator in J = 3 large A'^ QCD 



In order to test the proposed conjecture in section 3.2.2 we need a definition of the Wilson loop 
operator that does not suffer from perimeter divergence. We know that numerical computation 
of the string tension are performed by using correlators of smeared Polyakov loops. Therefore, 
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we define smeared rectangular Wilson loops of size nxm 

parameter / is the APE smearing factor and k is the number of smearing iterations, k should be 
proportional to the square of the perimeter for dimensional reasons since the smeared propagator 



is of the form [ 27 1 



V(^)=/(^)(5.v-^)+^ (4.5) 
f{q)=e'^^' (4.6) 



4.6.1 Test of the critical behavior of Wilson loops 



Wilson loops show critical behavior even without smearing []16|]. The distribution of the eigen 



values of the Wilson loop operator show good agreement with the Durhuus-Olesen distribution as 



given by (3.2) and (p^ 



It is necessary to study the continuum limit of this critical behavior and for this purpose one 
has to show that the eigenvalue distribution of W{n,m;f;k = ("+"') ;^;jv) in the Oc phase at a fixed 
n,m and / undergoes a transition from having no gap at small b to having a gap at large b as N °°. 
Furthermore, the critical bc{n,m;f;N) should scale properly as N ^ °°, nm —>■ oo such that 

lim ^= = — -— (4.7) 

has a finite limit. is the critical size of the Wilson loop and it will depend upon /. 

We fix the lattice size L and size of the color group, N. We then pick one value of / and pick 
a square Wilson loop, n = m. We study the eigenvalue distribution, e'^ of the Wilson loop operator 
for a range of b such that distribution goes through the transition. This is illustrated in Fig.^ where 
the eigenvalue closest to n is plotted as a function of for a 6 x 6 loop at N = 37 on a 8^ lattice. 



The non-zero gap is estimated [^] by matching the mean and variance to the Tracy-Widom [ |28| ] 
distribution for the universal distribution of the largest eigenvalue in the Gaussian ensemble. The 
explicit equation for the gap is 



i-i 

71 



(0^) + 1.964OO484J(02)_(0^)2 



(4.8) 



where On is the eigenvalue closest to n. 

4.6.2 Extracting the critical coupling and related parameters 

We perform a numerical calculation of the expectation value of 

ON{r,b) = (detiei+e^m)'^ (4.9) 

where W G SU (N) is a n x n Wilson loop at a lattice coupling b. We expect this observable to 
exhibit critical behavior atb = bc{n) and r = as ^ oo. If the double scaling limit is in the same 
universality class as the one in 2d then we expect 

i,m^-(;,,«)o„ (.= (^)'-|^.;, = ;,(„) + -^] =a?,„) (4.iO) 
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Figure 3: Distribution of the eigenvalue closest to n 



where ■A'{b,N) is a normalization factor and , a) is the Generalized Airy integral as given by 

A test of the above conjecture proceeds by first obtaining an estimate for bc{n), a\{n) and 
a2{n). Since det W = 1, it follows that ON{r,b) is an even function of r. It is also evident from 
that (^((^ , a) is an even function of . Let 



ON{r,b)=Co{b,N)+Ci{b,Ny + C2{b,N)r^ + ' 
be the Talyor's series for Ot^{r,b). Consider 

Co{b,N)C2{b,N) 



Q.{b,N) 



Cj{b,N) 



(4.11) 



(4.12) 



It is clear that this quantity will be the same for ON{r,b) and ^{b,N)ON(r,b). It is also clear that 

this quantity remains the same if we replace r by (3^) ■* ^^^tj ^rid view the Taylor's series as an 
expansion in ^. Therefore, the value of this quantity is well defined at the critical point, b = bc{n) 
and is given by 



0.364739936 



(4.13) 



The first equality in ( |4. 13| ) is obtained by evaluating the same quantity starting from ([3.9|). and we 
have used 



ir 

2 



2k+l 
4 



(4.14) 



Therefore, we obtain an estimate of Ci{b,N); / = 0, 1,2, using Montecarlo simulations and thereby 
obtain an estimate of Q.{b,N). We then use ( [^.13 ) to obtain an estimate of bc{n) at a fixed A'^ and 
extrapolate it to A/^ ^ 00. 
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The parameter a2{n) is defined via 

a2{nWN 

Tlierefore, we use tlie following relation 

dn{b,N). _ I da. _ r2(i) /r^^ 



b = bc{n) + —^^-^ (4.15) 



4f - 1 = 0.0464609668 (4. 16) 



da a2{n)VN db ^'='^^"' eVln \ 16k^ 

to obtain a2{n) at a fixed A'^. Since this is a sub-leading quantity, errors are larger in this quantity 
than in bc{n). 

Upon substitution of 

4 \3 ^ 



3Ny ai{n) 



in ( |4. llD , we conclude that 



4 1 Ciibcin),N) n 



3N^a\{n)Co{bc{n),N) ^T^{ 



(4.17) 



0.16899456 (4.18) 



and we use this relation to obtain an estimate of a\{n). The results as a function of N can be 
extrapolated to get the value SitN = °°. 

As an example of the above procedure, one finds be = 0.8095(4), a2 = 2.76(27) and a\ = 
0.8891(12) for a 4 X 4 Wilson loop at = 47 on a 8^ lattice with / = 0.03. The resulting function 



on the lattice as defined by ( [4.10| ) matches quite well with the Generalized Airy integral. Due 
to the arbitrary normalization that is involved in the matching, one possible way of checking the 
agreement is to look at the ratio, ^j^'^j for several values of a as a function of t, . Fig. ^ shows 
such a comparison. The approach to the large Umit of be, a2 and a\ are shown in Figs. ^-0. The 
agreement with the Generalized Airy integral gets better as one gets closer to the large limit. 
We end the discussion of large A'^ 3d QCD with a pictorial summary shown in Fig. |8[ 



5. Large A^ Ad QCD 

Large A^ Ad QCD was reviewed in Lattice 2005 [p9|]. We begin by summarizing the current 
status and focus on topics that were not covered in the previous review. 

There is a transition in the plaquette operator. This occurs at = 0.36 for all if' > 4^ and the 
transition is first order [ ^ , Both Oc phase and Ic phase have a continuum limit |32]. The 
critical size scales according to 

Le{b) = {0.250 ±0.025) ^^^J (5.1) 
This transition is the deconfinement transition studied on the lattice by taking the large A'^ limit 



using A'^ = 2,3,4, 6, 8 [33]. The deconfinement transition is first order and the latent heat has been 
measured through the jump in the internal energy. The latent heat is found to be de ^ 0.26A^^esB 
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Figure 4: The Generalized Airy integral is compared to lattice data for a 4 x 4 Wilson loop at A' = 47 on 8^ 
lattice. 
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Figure 7: Plot of ai {N) as a function of for 4X4 loop on 8^ lattice. 
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Continuum limit 



Oc -phase: Plaquette distribution 
opens up a gap around 9-7C. ^ 
QCD in the confined phase. Np° 
No volume dependence ^^d^ 




3c-phase: U(l) in all three directions are broken. 
QCD in a small box at high temperatures. 



Tadpole improved 't Hooft Coupling b 



Figure 8: Summary of large N QCD 'md = 3 onl? lattice 



where esB is the blackbody energy density per massless vector particle [p4|]. Critical sizes associ- 
ated with the lc-2c, 2c-3c and 3c-4c transition have not been determined yet. 



Transitions in smeared Wilson loops were first studied in Ad pq ] before starting the careful 
investigation of the double scaling limit in 3d. The transition fits the Durhuus-Olesen behavior as 
shown in Fig. ^. A careful determination of the critical area still needs to be performed using the 
double scaling limit. 

It is interesting to look at the fermionic sector of large A'^ Ad QCD. Chiral symmetry is broken 



in the Oc phase and the chiral condensate is found to be | |35| ] 

1 



{^r^rr\2G&W) k, (174MeV)' 



(5.2) 



Assuming = 3 is large enough, we get {^^)'^^ [IGdW) k, (251MeV)3 for SU(3). m\ oc niq as 



expected and [36| 



N 



;71MeV. 



(5.3) 



This translates to /;;: = 123 MeV for SU(3). This is the first instance we know of where Montecarlo 
simulations have indicated large 1 /N corrections. Pseudoscalar masses as well as vector meson 
masses were recently computed [ |39| ] for = 2,3,4,6 and extrapolated to the large A'^ limit. It 
would be interesting to study current correlators and compute vector meson masses directly in the 
large A'^ limit using quenched momentum techniques. It would also be interesting to study the 
correlations of Dirac eigenvalues as a function of force-fed momentum and their relation to fji 
using recent ideas from random matrix theory pO|]. 

The behav- 



Chiral symmetry is restored in the Ic phase and it is a first order transition [37, 3^ 
ior of the Dirac spectrum in the 3c and 4c phases would shed some useful insight into dimensionally 
reduced theories. 
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Ql , \ < \ < 1 , 1 I '^ l IN 1 

0.2 0.4 0.6 0.8 1 

9/Jl 

Figure 9: Fit of the lattice data to the Durhuus-Olesen distributions for four different sizes of wilson loops, 
namely, l„tc = 0.740,0.660,0.560,0.503. The associated areas (f in the 2d notation) that describe the con- 
tinuous curves are given by t = 8.06, 4.60, 2.82, 2.30 respectively. 



5.1 Twisted Eguchi-Kawai formalism 

Twisted Eguchi-Kawai was originally proposed as a way to preserve the U^{\) symmetries on 
a 1^ lattice all the way to the continuum limit pi]]. The basic idea is to modify the gauge action on 
the 1^ lattice to 

4 

5tek = -bN £ Tr (l - e-''^">^^U^UyUlul) . (5.4) 

The twist factors, n^y = —riy^ are integers. If one chooses, n^y = L; jJ. > V, and N = L}, then 
the theory behaves as if it is on a lattice. This theory has a symmetry like the usual Eguchi- 
Kawai model. But this symmetry is not broken in the strong coupling limit or weak coupling limit 
for the above choice of twists. 

There has been a recent revival of the twisted Eguchi-Kawai formalism of large Ad QCD 
with the aim of numerically studying the status of the symmetry as a function of b from b = Oio 
b = oo. A surprising new result is that twisted Eguchi-Kawai models on l'^ lattice seem to break the 
Z'^{N) symmetries for large enough N for a certain range of couplings [p^]. One sees a cascade of 
transitions where one goes from Z^^' ^ — > Z^ — > Z^ ^ Z^ ^ Z^ — > Z^^'^ . Z'^''^ and Z^^'^ corresponds 
to the fully symmetric phase at b = and b = oo respectively. The cascade of transitions do not 
occur for small A'^ (A'^ < 8 1 ) and begin to occur as bifurcations for A/^ > 8 1 . These transitions seem to 
be strongly first order. For large enough coupling at a fixed A'^, all the Z"* (A'^) are most likely restored. 
It is hard to see the restoration numerically starting from the fully broken phase and this is also the 
reason why one cannot confirm if there is a reverse cascade that takes it from Z^ Z^'^ . The 
above result has been confirmed independently by Ishikawa |43] and Okawa [44]. The dependence 
on the critical coupling for these transitions as a function of L has not been studied. The possibility 
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to define a continuum limit of large N QCD using twisted Eguchi-Kawai model will depend on 
how the critical coupling scales with L. Also, other twists, perhaps with prime numbers for A'^ may 
show a different picture. 

6. Theta parameter in large Ad QCD 

In the instanton dilute gas approximation, the free energy as a function of the d parameter is 



given by [45] 



f{e)=x{i-co^e) (6.1) 

where % is the topological susceptibility which is expected to go to zero exponentially in A^. Large 
A'^ arguments suggest that Q /N is the parameter to keep fixed as A'^ ^ oo. Since f{Q) should be a 
periodic function in 6, f{Q) cannot be an analytical function of 6: 

f{e) = ^rmn{d + 2nkf (6.2) 

2 k 

Instanton and large A^ make qualitatively different predictions for the moments of the topolog- 
ical charge. Instanton arguments say that the topological susceptibility goes down exponentially 
with A^. Large A^ arguments give a finite topological susceptibility in the large A'^ limit. The ratio of 
the fourth moment to the second moment would be unity for instantons and would be zero by large 
A'^ argument. Large A^ predicts 

lim < 2 > |;r-e / lim < 2 > 


Lattice computations [BtI] of the topological susceptibility for A^ = 2,3,4,6,8 show that the 
large A'^ limit is finite and is given by = 0.390(14). A high statistics computation of the topo- 
logical charge for SU(3) [48] and a field theoretical approach [^] show that the ratio of the fourth 
moment to the second moment is significantly smaller than unity. Both results favor the large A'^ ar- 



gument. Gauge theories have also been studied on the lattice by an expansion around 6=0 [ pO , 
These results also are in agreement with the large A'^ predications. A direct measurement of the non- 
analyticity at = tt is difficult. 

7. Principal chiral models in J = 2 

Two dimensional principal chiral models for large A'^ are similar to four dimensional large A'^ 



gauge theories. The principal chiral model is defined through the action [ ]52[ ] 



S = -NbYj^vp {x) {u\x + (i) + U\x-fL))] (7.1) 

This model has a global SU (N) xSU (N) symmetry under which 

Uix)^R'^U{x)L; /?,LgSU(N) (7.2) 

and it undergoes a second order phase transition at be = 0.3057(3). The theory is in the continuum 
phase for b > be and the continuum Umit is reached by taking b ^oo. 
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Define the operator 



L{n,m;b) =U{x)U^{x + njX+mv)', n<m. 



(V.3) 



It can be used to define the correlation function 



G{n,m;b) = — {TrL{n,m;b)) 



The correlation length is defined as 



4sin2f 



G{0,l;b) 



(7.4) 



(7.5) 



where G{p\,p2',b) is the lattice Fourier transform of G{n,m;b). The correlation length diverges at 
b = be and also as b ^oo. 

Like in the case of large gauge theories, it is useful to consider the eigenvalues of L{n,m;b). 
These eigenvalues are invariant under the global transformations. Let p{Q;n,m;b) define the dis- 
tribution of eigenvalues. L(0, \ ;b) appears in the action and is analogous to the plaquette operator 
in large N gauge theories. p(0;O, \;b) does not have a gap for b <bc and has a gap for b > be [^]. 
The universal behavior has not yet been analyzed. 

/ — ^ 1 

L{n,m;b) is analogous to a Wilson loop operator in large n gauge theories with r = "^^^^ 
being the physical length. We expect p{6;r) in the continuum limit to show critical behavior such 
that it has a gap for r < rc and it does not have a gap for r > r^. An initial investigation [53] 



of the gap as defined in section [4.6. 1| is plotted as a function of r in Fig. [10| and it suggests the 
expected picture. But, a closer look indicates a drift in the critical value of r as one gets closer to 
the continuum limit. This might be an effect of not using smeared U{x). A proper investigation will 
have to use smeared U {x) and the critical size will have be studied as a function of the smearing 
factor. 



8. Large gauge theories with adjoint matter 

Large gauge theories on a <i dimensional torus with p adjoint matter fields can be viewed as a 
{d + p) dimensional large gauge theory in the pc phase where the length of the broken directions 



are taken to be zero[54, 55]. Let the masses of the p scalar fields be the same and let the lengths of 
the periodic directions of the d dimensional torus be free parameters of the theory. The Polyakov 
loops associated with the gauge fields on the d dimensional torus serve as order parameters. Such 
theories 'md = \ and d = 2 can be analyzed pseudo-analytically (with a little bit of numerical help) 
for some region of the parameter space. 

The single Polyakov loop m d = \ breaks if p > \. Pseudo-analytical analysis shows that 
there is a line in the two dimensional coupling constant plane (mass and length of the circle) that 
separates the broken phase from the unbroken phase. This is consistent with the existence of 2c to 
3c phase transition in J = 3 and a 3c to 4c transition m d = 4 for large A'^ gauge theories onu d 
dimensional torus. 

There are two Polyakov loops and p>\. It has only been possible to perform a large mass 
(of the adjoint scalar field) analysis. The large mass analysis predicts three phases: (i) Both loops 
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Figure 10: Gap as a function of r in the distribution of p(0;r) 



are not broken; (ii) One of the loops is broken; (iii) Both loops are broken. This result is again 
consistent with the existence of Ic ^ 2c — > 3c cascading phase transition in J = 3 and a 2c — > 
3c — > 4c cascading transition in J = 4 for large N gauge theories on a J dimensional torus. 



9. Gregory-Laflamme transition 

Consider a d dimensional theory of gravity with no matter fields and consider the background 
spacetime to be x T". Assume the torus to be of the same length, L, in all directions. A 

p-brane solution is a Schwarzschild black brane that is independent of p of the n directions on the 
torus, p = is a black hole. The entropy of the black brane defines a temperature, Th, and t = ThL 
defines the length of the torus in terms of the black brane temperature. There exists a.t = ?gl(p) 
such that p-brane decays into a (p — 1) -brane as t increases through ?gl(p) this referred to as 
the Gregory-Laflamme transition [^6|]. ?gl(p) < ^GL(p-i) ^rid there exists a fc(p-i) such that ?gl(p) < 
^C(p-i) < ^GL(p-i)- The free energy for a (p — 1) brane is favored to a p brane as t increases through 
^c(p-i)- This cascade of transitions is like the Oc ^ Ic ^ • • • ^ Jc cascade observed in large 



gauge theories \p%. There is a relation between these two transitions - Fermions can be discarded 
in super Yang-Mills at high temperatures since they obey anti-periodic boundary conditions and 
the theory reduces to Yang-Mills with adjoint scalars. 



10. Other related topics 

There are several other recent developments in the area of a large N gauge theories that were 
not presented due to time constraints. 

Several papers considered the case of fermionic matter in two-index representations, adding 
order degrees of freedom and consequently changing the large N dynamics of the pure gauge 
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field 1 58, 59, S^]- In particular, the extra repulsion between the eigenvalues of Polyakov loops the 
matters fields perturbatively generate can delay or remove the bulk transitions of the pure gauge 
system. In one case there is an argument for the absence of all bulk transitions, indicating volume 



independence down to zero size in the continuum []59|]. We do not know of any numerical work 
testing this prediction. 

Another topic is the addition of chemical potential for the quark fields. Here again, fermions 
will play a dynamical part and we are not aware of numerical work pertaining to the large limit. 
Some discussion of the physical implications of the chemical potential in large A'^ QCD can be 



found in [61] and [52] 
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